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A theoretical study ol the transport properties ol zigzag and armchair graphene nanoribbons 
with a magnetic barrier on top is presented. The magnetic barrier modifies the energy spectrum ol 
the nanoribbons locally, which results in an energy shift of the conductance steps towards higher 
energies. The magnetic barrier also induces Fabry — Perot type oscillations, provided the edges of the 
barrier are sufficiently sharp. The lowest propagating state present in zigzag and metallic armchair 
nanoribbons prevent confinement of the charge carriers by the magnetic barrier. Disordered edges 
in nanoribbons tend to localize the lowest propagating state, which get delocalized in the magnetic 
barrier region. Thus, in sharp contrast to the case of two-dimensional graphene, the charge carriers 
in graphene nanoribbons cannot be confined by magnetic barriers. We also present a novel method 
based on the Green's function technique for the calculation of the magnetosubband structure, Bloch 
states and magnetoconductance of the graphene nanoribbons in a perpendicular magnetic field. 
Utilization of this method greatly facilitates the conductance calculations, because, in contrast to 
excising methods, the present method does not require self-consistent calculations for the surface 
Green's function. 

PACS numbers: 73.23.Ad, 75.70.Cn, 73.63.Bd 



I. INTRODUCTION 

The single planar sheet with carbon atoms densely 
packed in a honeycomb structure forms the so-called 
graphene, which demonstrates a variety of unique elec- 
tronic transport properties and has the potential ap- 
plications in the future nanoclectronicsi. Theoretical 
studies have indicated that the special lattice structure 
of the graphene results in nearly linear dispersion rela- 
tions around the K points (Dirac points) of the Bril- 
louin zone^. This unique band structure is responsible for 
the distinct electronic properties of the graphene. Near 
the Dirac point, electrons manifest themselves the mass- 
less chiral fermions and can be described by the Dirac 
equatior>2ii£. The electronic transport behaviors of the 
two-dimensional graphene subjected to an electrostatic 
potential^ or a magnetic barrier (MB)£ were studied on 
the basis of the Dirac equation, which indicate that the 
Dirac fermions can be transmitted perfectly through a 
classically forbidden region while confined effectively by 
the magnetic barrier. Moreover, the anomalous integer 
and fractional quantum Hall effects in two-dimensional 
graphene have been studied experimentally and theoret- 
ically by various groupsi^^. 

The rollcd-up graphene is known as the single-wall 
carbon nanotube whose electronic properties have been 
studied extensively in the past decades. The quantized 
conductance and Fabry — Perot interference pattern were 
observed experimentally and interpreted by various the- 
oretical approaches^. The other interesting effects in- 
cluding Coulomb blockade^ and Kondo effects^, and 
the electronic transport in ballistic^ and disordered 



nanotubesi^ were studied. Another related carbon-based 
structure is the graphene nanoribbon (GNR), referred to 
the quasi-one dimensional graphene with a finite width 
W. Recent development of the experimental technique 
enable one to fabricate very narrow GNRs with ultra- 
smooth edges of the width W < 10 nmi&. The elec- 
trons propagate in such narrow systems very differently 
compared with the two-dimensional graphene where the 
edges are totally irrelevant. In graphene ribbons, the 
transport properties are strongly influenced by their 
edges along the transport direction which are distin- 
guished into two types: zigzag and armchair. For arm- 
chair case, it is particularly interesting that the graphene 
ribbons may be metallic or semiconducting depending on 
their widths. There is a lot of theoretical effort devoted 
to the studies of the quantum transport in graphene 
ribbons. The conductance quantization in mesoscopic 
graphene^ and coherent transport in graphene nanocon- 
strictions with or without defects^ were reported re- 
cently. 

The purpose of the present paper is twofold. First, we 
explore a possibility to control electron conductance of 
graphene nanoribbons with the help of magnetic barri- 
ers. MBs in the conventional quantum wires (QWRs) 
have been the subject of theoretical and experimental 
studies, which are driven by the MB's potential abil- 
ity of parametric spin filtering. The pioneering theo- 
retical research by Peeters et al»i& indicated that the 
magnetic barrier possesses the wave- vector filtering prop- 
erties in QWRs and further work in graphene was also 
suggested^. Furthermore, recent theoretical studies have 
revealed further rich phenomenology of magnetic barri- 
ers in quantum wires, such as Fano-type resonances^ and 



spin filterin g 20 ! 21 . In two-dimensional graphene, theoret- 
ical work has shown the strong effects of the magnetic 
barrier on the direction-dependent transmission^. Our 
studies will focus on the magnetic barrier effects on the 
quasi-one-dimensional GNRs. 

Second, we present a detailed description of a novel 
method based on the Green's function technique for 
the calculation of the magnetosubband structure, 
Bloch states and magnetoconductance of the graphene 
nanoribbons in a perpendicular magnetic field. Note 
that magnetoconductance calculations for the graphene 
nanoribbons based on the Green's function technique 
has been reported previousl y 17 ' 22 . However, a distinct 
feature of the present method is a novel approach 
to calculation of the surface Green's function T for 
semi-infinite nanoribbons. In contrast to the Green's 
functions for finite structures that can be easily calcu- 
lated by adding slice by slice in a recursive way with 
the help of the Dyson's equation, the calculation of 
the surface Green's function of a semi-infinite structure 
represents a non-trivial problem. Such calculations are 
typically done self-consistently which makes conductance 
calculations very time-consuming. In the present paper 
we present a different method of computing T which 
does not require self- consistent calculations. Instead, 
the surface Green's function is expressed via the Bloch 
states of the graphene nanoribbons which in turn are 
simply obtained as solutions of the eigenequation of 
the dimension 2N x 2N (with N being the width of 
the nanoribbon). Utilization of this method greatly 
facilitates the conductance calculations, making the 
present method far more efficient in comparison to the 
existing ones. Programming codes for calculation of the 
magnetosubband structure, the surface Greens function 
and the magnetoconductance based on the developed 
method are freely available in the AIP EPAPS electronic 
depository^ 

This paper is organized as follows: In Sec. II we sketch 
the geometry of the devices and briefly introduce the 
model for the conductance for our calculations. In Sec. 
Ill, we describe the tight-binding model for the graphene, 
theory of the Green's function method, as well as the for- 
malism for the computation of surface Green's functions. 
This is followed by the presentation and discussion of the 
numerical results in Sec. IV. Summary and conclusion 
constitute the Sec. V. 



II. FORMULATION OF THE PROBLEM 

The geometries under consideration for graphene 
nanoribbons with zigzag and armchair edges are illus- 
trated in Fig. 1(b) and (c), respectively, where the left 
and right leads are made of semi-infinite graphene. The 
nanoribbons are subjected to a magnetic barrier whose 
shapes may be rectangular or smooth as shown in Fig. 
1(a) with zero magnetic field in leads. (Note however 
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FIG. 1: Schematic geometry of the structure under considera- 
tion for the case of (b) zigzag and (c) armchair graphene. The 
current through the central part of the device is injected and 
collected in semi-infinite ideal leads representing graphene 
nanoribbons of width N. Unit cells of the graphene nanorib- 
bons are marked by blue dashed rectangles (see Fig. 0. (a) 
A ferromagnetic film deposited on the top of the graphene 
nanoribbons gives rise to an inhomogeneous magnetic field. 



that the theory presented in the next section is not re- 
stricted to the case of zero field in the leads). The mag- 
netic barrier represents a strongly localized magnetic field 
that is oriented perpendicular to the surface of the rib- 
bon. Magnetic barriers with amplitudes up to IT have 
been realized experimentally by ferromagnetic films on 
top of a graphene shee t 24 ' 25 : magnetizing the ferromag- 
netic film in the transport direction results in a magnetic 
fringe field with a perpendicular component localized at 
the edge of the film that extends along the transverse 
direction. Alternatively, magnetic barrier formation has 
been demonstrated by placing two-dimensional electron 
gases with a step in an external magnetic field j2£ an ap- 
proach which conceptually allows much larger barrier am- 
plitudes. Both concepts should be in principle adaptable 
to graphene nanoribbons. 

We model the leads and the device in the middle by 
the standard tight-binding Hamiltonian on the honey- 
comb lattice, see below, Eq. ([3]). The conductance G 
can be calculated using the Landauer-Biittiker formal- 
ism which gives the conductance of the system in terms 
of the electron transmission coefficient T. expressed as 
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FIG. 2: Geometry of (a) zigzag and (b) armchair graphene 
ribbons. The ribbons are periodic in the i-direction (slices 
in the ^-directions are labeled by index i). The figure shows 
nanoribbons with N — 6 sites in the transverse direction. 
Unit cells of the graphene nanoribbons are marked by dashed 
rectangles. 



The standard tight-binding Hamiltonian has the form 

H = V r a+a r - t r . r+ Aa+a r+/ \, (3) 



r.A 



where V r describes the electrostatic potential on the site 
r = i,j and summation in the second term is performed 
over all available nearest neighbors with i r , r +A being the 
nearest-neighbor hopping integral. In the absence of a 
magnetic field the nearest-neighbor hopping integral is 
tr,r+A = to ~ 2.7 eV. In the presence of an external per- 
pendicular magnetic field B the hopping integral acquires 
the Peierls phase factor, t r _ r+/ \ = to exp(i0 riI ._|_A),where 
&r.r+A = ^4>r.r+A/4>o j with 4>r.r+A being the line inte- 
gral of the vector potential A from site r to a neighboring 
site r + A, 



r+A 



A-dl 



(4) 



and <^o = h/e is the flux quantum (in our calculations 
we use the Landau gauge, A = (— By,0)). (In calcula- 
tion of hopping integral ((3]) we use the carbon-carbon 
bond length a = 0.142 nm, see Appendix). Note that 
the Hamiltonian operator H is convenient to write down 
in the form, 



where T{E) is the total transmission coefficient, Jfd(E— 
Ep) is the Fermi-Dirac distribution function and Ep is 
the Fermi energy. 

We calculate the transmission amplitudes of electrons 
injected to the systems using the recursive Green's func- 
tion method which is described in the next section. 



III. THEORY 
A. Basics 

We define the Bloch states in the infinite periodic 
graphene ribbons, 

I i>) = E ^ <3 1 °) • ^ = ^ kx * > ( 2 ) 

where afj (ciij) is a standard creation (annihilation) op- 
erator on the site ipij is the amplitude of the wave 
function on the site Xi is the coordinate of the i-th 
slice, k is the Bloch wave vector in the direction of the 
translational invariance x, and the summation runs over 
all sites of the graphene lattice (see Fig. [J). Note that 
this form of the wave function does not distinguish be- 
tween sublattices A and B of the graphene lattice. An 
explicit distinction between these sublattices is not nec- 
essary when using the Green's function technique, where, 
instead, it is more convenient to define the wave function 
on slices of the lattice (see Sec. IIIIB[) . 



ff = E [hi] + U, (5) 

i 

where hi describes the Hamiltonian of the z-th slice, and 
U describes hopping between all neighboring slices (ex- 
plicit forms of hi and U can be easily obtained from Eq. 
©)■ 

The Green's function of the operator H is defined in a 
standard wa y 27 i 28 , 

(E-H + ie)G = I (6) 
where / is the unitary operator. 



B. Bloch states and velocities in the graphene 
nanoribbons. 

We continue by describing a method for calculation 
of the Bloch states and their group velocities in the 
zigzag and armchair graphene nanoribbons in the pres- 
ence of a perpendicular magnetic field. The method is 
based on the technique developed for calculation of the 
band structure of a mesoscopic antidot lattice in confined 
geometries 2 ^ and has been used for calculation of the 
Bloch states in photonic structures^ and in the interact- 
ing quantum wires in the integer quantum Hall regime^. 

Consider an infinite ideal graphene ribbon with N sites 
in the transverse j-direction, Fig. [2] A unit cell of the 
structure consists of M slices, where M = 2 for the zigzag 
graphene and M = 4 for the armchair graphene. 
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The Hamiltonian of an ideal infinitely long graphcnc 
ribbon can be written in the form 



Combining Eqs. 
eigenequation, 



T2l) and (Tl"3l). we arrive at the 



H = H, 



cell 



(7) 



where the operators H ce \\ and H out describe respectively 
the unit under consideration (1 < i < M), and the out- 
side region including all other slices — oo < i < and 
AI + 1 < i < oo, and [/ is the hopping operator between 
the cell and slices i = and i = M + 1 (an explicit form 
for these operators can be easily obtained from Eq. ([3])). 
We write the total wave function, Eq. ([2|), in the form 



IV-') = l^ccll) + |Vw), 



(8) 



where iV'ccii) and li/jout) are respectively wave functions 
in the cell and in the outside region. Substituting Eqs. 
(HI,© into the Schrodinger equation H\ip) = E\tjj) and 
using the definition of the Green's function, Eq. (JS]), we 
obtain IV'coii) = G co iiJ7|'!/'out), where G C eii is the Green's 
function of the operator H ce u. Taking the matrix ele- 
ments of the wave functions in the real space represen- 
tations, = (O a ij|V') f° r the nrs t (i = 1) and the 
last (i = M) slices of the unit cell, this equation can be 
written in the matrix form, 



V'i = G'ccn^i.oV'o + Gila U? fi il> M +i 
V'm = G™JU lfi ip + G^fU+ipM+i, 



(9) 



where tpi is the vector column describing the wave func- 
tion for the slice i, 



A = (Vv; ■ ■ ■ ! V'i,A r ) J 



(10) 



and C/i.o and G^'cii denote the matrixes with the matrix 
elements 



(Ui,o) jjt = (0a hj \U\a+ r 0) 



(11) 



(Gecn) 



j j 



(Oa^lGcenla+,,0). 



Explicit expressions for the matrix elements of the matrix 
U are given in the Appendix. In the derivation of Eq. 
© we used Um.m+i = C^o.i (because of the periodicity 
of the ribbons) and Uo.i — U^ ('+' stands for Hermitian 
conjugate). 

It is convenient to rewrite Eq. ([9]) in a compact form 



1pM+l 



= T 2 



01 
00 



, where 



(12) 



1 _ [-Gsrm* 



T 2 



o G%£ Ult0 



with / being the unitary matrix. The wave function of 
the periodic structure has the Bloch form, 



VVn+M 



AkM 



(14) 



determining the set of Bloch eigenvectors k a and eigen- 
functions tp a , 1 < a < N. It should be stressed that this 
eigenequation provides a set of the Bloch states {k a } for 
a fixed energy E, which includes both propagating and 
evanescent states. The latter can be easily identified by 
a non-zero imaginary part. 

In order to separate right- and left-propagating states, 
fc+ and k~ , we compute the group velocities of the Bloch 
states v a = ■§§-, whose signs determine the direction of 
propagation ('+' stands for the right-propagating and '- 
' for the left propagating states). The group velocities 
can be computed directly by numerical differentiation of 
the dispersion relation. This is however not an efficient 
approach because for each energy the eigensolver gives 
cigenstatcs a in different order. We instead derive below 
a simple formula which gives the group velocities of the 
Bloch states based on the eigenfunctions of Eq. (fT4|) . 

Consider a unit cell of an infinite graphene nanoribbon 
consisting of M slices. The wavefunction of the a-th 
Bloch state can be conveniently rewritten in the form 

V>) = S»=i iV'i)) where \ipi) is the wave function for the 
i-th slice, 



<Pi) 



(15) 



[To simplify our notations we have dropped the Bloch 
index a ]. Starting from the Schrodinger equation and 
calculating the matrix element of the Hamiltonian of 
the unit cell, we obtain for each slice i, (^i\H\tp) = 
E(ipi\ip) = E\ipi\ 2 . Performing summation over all slices 
of the unit cell and using a definition of the group veloc- 
ity, we obtain 



dE 
~dk 



1 



AI 



d 



M ^ dk 



(Am: 



m\ 



(16) 



where the summation is performed over all slices i of the 
unit cell, and 

<Pi = (<Pi,i;-;<Pi,N) T (17) 

is a vector composed of the matrix elements (fij = 
(Oai t j\tp) .(Note that according to Eq. (TiT)|) and iflB"]). 
vectors (fi can be obtained from ipi via the relation 
V'i = e lkXi (fi). Representing the Hamiltonian of the unit 
cell in the form ([5]), the matrix elements (ipi\H\ip) can be 
easily evaluated, which gives 



. AI 

-y 



■•Pi 



(xi - Xi-i)U it i-iipi-ie 



ik{xi 



-(xi+i - Xi)Ui li+ i<pi + ie lk{x * +1 Xl) 



i) 



(18) 



where the matrixes L^y are defined by Eq. (jlip [explicit 
expressions for these matrix elements are given in the 



(13) Appendix]. 
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C. Surface Greens function V. 

Here, we describe an efficient method for calculation 
of the surface Green's function T in the magnetic field^ 
Note that most of the methods for calculation of the 
Green's function reported to date require searching for 
a self-consistent solution for T which makes these calcu- 
lation very time consumin g 17 ' 22 . In contrast, our method 
does not require self-consistent calculations, and the sur- 
face Greens function is simply given by multiplication of 
matrixes composed of the Bloch states of the graphcnc 
lattice (see below, Eqs. (f2"Tj) , (j2"2"]) ). The calculations de- 
scribed in this section are based on the method developed 
in Refi^ for periodic photonic crystals which is adapted 
here for the case of the graphcnc nanoribbons. 

Consider a semi-infinite periodic ideal graphcnc ribbon 
extended to the right in the region — m < i < oo. Sup- 
pose that an excitation \s) is applied to its surface slice 
i = —in. Introducing the Green's function of the semi- 
infinite ribbon, G r ib, one can write down the response to 
the excitation \s) in a standard form 2 ^. 

\ip) = G rib |s}, (19) 

where \tp) is the wave function that has to satisfy the 
Bloch condition Consider a unit cell of a graphene 
lattice, 1 < i < M, (M = 2 and 4 for the zigzag and 
armchair lattices, see Fig. Applying Dyson's equation 
between the slices and 1 we obtain 

G 1 I i- m =T r U li oG° I £ m , (20) 

where T r = G ri ' b is the right surface Green's function 
(i.e. the surface function of the scmiinfinitc ribbon open 
to the right), and the definition of the matrixes U and G 
in the real space representation are given by Eq. (jTTJ)). 
Evaluating the matrix elements (Oai j of Eq. ([!§)) and 
making use of Eq. (|20[) . we obtain for an each Bloch state 
a, ipf — T r Ui.oipQ. The latter equations can be used for 
determination of r r , 

r r u lt0 = *i* \ (21) 

where and are the square matrixes composed of 
the matrix-columns tpf and ipQ , (1 < a < N), Eq. (TT4")) . 
i.e. = {^i,...,ipi); *o = C0o,-,^oO- The expression 
for the left surface Greens function Ti (i.e. the surface 
function of the semiinfinite ribbon open to the right) is 
derived in a similar fashion, 

TtU+ Q = ^ M ^ +1 , (22) 

where the matrixes $ m and ^m+i are defined in a sim- 
ilar way as and <3/o above. Note that matrixes ^>m 
and 'I'm+i can be easily obtained from and using 
the relation (fl"2"|) . Note also that when the magnetic field 
is restricted to zero, the right and left surface Greens 
functions are identical, Ti = T r . 



D. Magnetoconductance of the graphene 
nanoribbons 

In order to calculate the transmission coefficient T(E) 
we divide the structure into three regions, two ideal semi- 
infinite leads of the width N extending in the regions 
i < and i > L respectively, and the central device 
region (where scattering occurs), see Fig. [T] We assume 
that the left and right leads are identical. The incoming, 
transmitted and reflected states in the leads, \ip a ), \ip a ) 
and have the Bloch form (J5J), 

N 

l^) = E e ^E^;<;l°> ( 23 ) 

i<0 3=1 

N 

\o = EEW fe ^-* i) E<; a i-i°> ( 24 ) 

i>L j=l 
N 

io=EE^ Ali E<i<ii°)> (25) 

i<0 j=l 

where tp a (rp a ) stands for the transmission (reflection) 
amplitude from the incoming Bloch state a to the trans- 
mitted (reflected) Bloch state (3, and we choose xq = 0. 
The transmission and reflection coefficients are expressed 
through the corresponding amplitudes and the Bloch 
velocities 2 ^ 

t = E^i 2 ;^ = E-W 2 > 

a,p a,f3 

where the summation runs over propagating states only. 
The transmission and reflection amplitudes can be calcu- 
lated from the equations^, 

$ 1 T = -G L > (U ,i* 1 K 1 -Tr 1 $ ) (26) 
$ i? = -G -°([/ ,i$i^i - rr^o) - $0 (27) 

where the matrixes T and R of the dimension N x N prop 
are composed of the transmission and reflection ampli- 
tudes (T)p a = tp a , (R)p a = r 0a ; (with N prop being 
the number of propagating modes in the leads); G L,a 
and G 0,0 are the Green's function matrixes with ma- 
trix elements defined according to Eq. (Till ; T; is the 
left surface Green's function, Eq. (|2"2"jk Uq,i is the hop- 
ping matrix between the left lead and the device region 
(jTTJ) ; K\ is the diagonal matrix with the matrix elements 
(K\) a fi = exp(ifc+ xi)5 a ,/3- The square matrixes $i and 
3>o describe the Bloch states on the slices 1 and of 
a ribbon unit cell (see Fig. \2§ and are composed of 
matrix-columns (j)" and </>q, (1 < a < N), Eq. (|17p . 

i.e. § 1 = (^,...,^);$o = (^ ) ...,^). 

Calculation of the Green's functions G L '° and G°'° is 
performed in a standard way 2 ^. We start from the Greens 
function of the first slice in the device region and, using 
the Dyson's equation, add recursively slice by slice un- 
til the last slice of this region is reached. Finally, we 
apply the Dyson's equation two more times adding the 
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left and right scmi-infinitc ribbons whose surface Green's 
functions are given by Eqs. (f2"Tj) . (f2"2"|) . 

Having calculated the transmission and reflection am- 
plitudes that give the wave functions on slices i = and 
i = L, we can easily restore the wave function inside the 
device region using the relation between the wave func- 
tions on slices i, V and i + 1, i' — 1 (we assume that i' > i) 

A+i = G^; +1 CA +M ^ + G^- 1 U+_ 1 ^, (28) 



inner 



where G^ el is the Green's function of the internal region 
only (extending from the slice i to the slice to). (Equation 
(f28|) is derived in a similar way as Eq. ((9])). Removing 
slice by slice from the inner region and repeatedly using 
Eq. ([28]) on each step, we restore the wave function in 
the entire region < i < L. 

The diagonal elements of the total Green's function for 
each slice i give the local density of states (LDOS) at the 

site p(i,j,E) = J . The LDOS can 



jj 



be used to calculate the local electron density at the site 



• J- 



n(i,j)= / dEp(i,j,E)f(E-E F ). 



(29) 



For quasi-one dimensional structures considered in this 
study it is convenient to introduce the local density of 
states integrated in the transverse direction, 



p(i,E)=J2p(hj,E). 



(30) 



IV. RESULTS AND DISCUSSION 

In this section we discuss the conductance properties 
of two-terminal GNRs with MBs using the formalism de- 
scribed above. GNRs with both zigzag and armchair 
edges are considered. The electronic properties of arm- 
chair GNRs depend strongly on its width W. The arm- 
chair GNRs are metallic when 2N + 1 is a multiple of 
3, and otherwise they are semiconducting. Metallic arm- 
chair GNRs behave similarly to zigzag GNRs regarding 
the effects discussed here, even though the origin of the 
first subband is different, and are not presented sepa- 
rately. 

Figure [3] shows the Fermi energy dependence of the 
conductance for the zigzag and armchair ribbons with 
N = 151 and 150, respectively, corresponding to a width 
of W w 32 nm. The rectangular magnetic barrier has a 
length of 120 nm. The smooth magnetic barrier has the 
standard shape realized in experiment a 25 i 33 and a full 
width at half maximum of 120 nm. For the case of the 
smooth barrier the central (device) region has a length 
of 360 nm. The shapes of the smooth and sharp barriers 
arc depicted schematically in the insets to Fig. [3J We 



present the conductance calculations for the maximum 
magnetic field strength in the barrier in the interval of 

— 8 T. While inhomogeneous fields up to 3.4 T have 
been achieved in the laboratory by using etch facets^ we 
consider such high fields in order to addre ss the regime 
when the magnetic length Ib = \Jh/eB (= 26 nm at 

1 T) is smaller than the ribbon width. Alternatively, this 
could have been achieved by increasing the ribbon width, 
which is however rather impractical from the computa- 
tional point of view. 

In the absence of MBs, the ballistic conductance of the 
GNRs is simply proportional to the number of subbands 
Nq at the Fermi energy at zero magnetic field^ see Fig. 
The conductance shows plateaus and increases as a 
function of Fermi energy, in analogy to the case of QWRs. 

Figures^a), (b) show the conductance of the semicon- 
ducting armchair GNR for the rectangular and smooth 
magnetic barriers. The dashed lines indicate the num- 
ber of propagating states Nb in the corresponding GNR 
in the homogeneous magnetic field whose amplitude is 
equal to the maximum field B in the barrier region. As 
the magnetic field increases the subbands depopulate and 
hence the corresponding number of available propagating 
states Nb decreases. Because the magnetic field provides 
an additional confinement in the ribbon, at a given Fermi 
energy the number of the magnetosubbands Nb is always 
smaller than No. Because of this, Nb represents the lim- 
iting factor for the conductance of the magnetic barrier 
structure such that No incoming states in the leads are 
redistributed among Nb available states in the magnetic 
barrier. This is clearly seen in Figs. [3] (a), (b) where the 
conductance of the structure at hand approximately fol- 
lows Nb ■ Note that the magnetic field reduces the energy 
gap in the vicinity of E = 0. Despite of this the conduc- 
tance of the magnetic barrier is always zero below the en- 
ergy threshold of E t h ~ 0.006in regardless of the strength 
of the magnetic barrier. This simply reflects the fact that 
propagating states are injected from the leads where the 
magnetic field is absent and the threshold propagation 
energy E t h ~ 0.006to is not affected by the strength of 
the barrier in the central region of the device. In ad- 
dition, transmission resonances are superimposed on the 
conductance plateaus. They are well pronounced for the 
rectangular barriers, but get heavily suppressed as the as 
magnetic barrier assumes the more realistic, soft shape. 
As the strength of the barrier increases, the resonances 
become more prominent. 

In the zigzag GNR with a MB, the conductance steps 
also move towards higher energies and follow Nb vs en- 
ergy as B increases, see Fig. EJc), (d). This, as in the 
case of the armchair GNRs, simply reflects the magnetic 
field induced shifts of the GNR modes in the barrier 
region. Around E = , an energy interval exists in 
which only the lowest propagating state contributes to 
the conductance. This state evolves from the disper- 
sionlcss edge state present in the zigzag GNRs at zero 
energy. The MB is thus able to reduce the number of 
current carrying states in certain energy intervals, e.g. 
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FIG. 3: The calculated conductance (full lines) and the number of occupied modes at the maximum magnetic field (dashed 
lines) as a function of the Fermi energy in the semiconducting armchair GNR (a,b) and the zigzag GNR (c,d) for MBs with 
different amplitudes 0T, 2T, 4T, and 8T. The inset in (a) represents the section of the energy dispersion at B — 8T which 
causes the change of the number of modes from 3 to 2 and back to 3 as the energy increases. 



between E « 0.03t and E w 0.037i for the MBs with 
a strengths of 8T. Note that for the zigzag GNR the 
conductance changes in steps of 2 x 2e 2 /h, whereas for 
the armchair GNR it changes in steps of 2e 2 /h. This 
reflects the difference in evolution of the subband struc- 
ture of corresponding homogeneous armchair and zigzag 
GNRs, where the number of states at the given energy 
depends on the wire width N and on whether the ribbon 
is metallic or insulating. (The conductance quantization 
for armchair and zigzag GNRs was discussed by Peres et 
al.)M 

In addition, as in the case of the armchair GNR, trans- 
mission resonances are observed for the rectangular MBs. 
These resonances are completely suppressed for the case 
of the smooth barriers. Both the frequency and the am- 
plitude of the oscillations become higher as the strength 
of the MBs is increased (Fig. [3]). Furthermore, the fre- 



quency decreases as the length of the MB is decreased 
(not shown). This behavior is similar to the conduc- 
tance resonances in quantum point contacts with abrupt 
openings^! and originates from multiple reflections at the 
edges of the MB along the transport direction. The mul- 
tiple reflections at the edges lead to the Fabry — Perot 
type oscillations, as can be seen in Fig. [4] for the case 
of a GNR with zigzag edges where the number of max- 
ima in the LDOS along the transport direction changes 
by one for successive resonances. Similar to the case of 
a smooth quantum point contacl^, a gradual change of 
the magnetic field reduces the reflection probabilities and 
suppresses the resonances, resulting in smaller oscillation 
amplitudes. 

The dependence of the conductance on the magnetic 
barrier suggests that complete confinement by magnetic 
barriers is not possible due to the presence of the low- 
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FIG. 4: The local density states (LDOS) with the rectan- 
gular MB of 8T integrated over the y direction (in units of 
t ( J 1 nm _1 ), as a function of energy and x in the zigzag GNR. 
The LDOS plot shows the localized states form inside the 
magnetic barrier. The position of the localized states in rela- 
tion to the corresponding conductance. 



est propagating state, in stark contrast to the case of 
two-dimensional graphene sheets^ To shed more light 
on the influence of the MB on the lowest propagating 
state, we study the local density of state (LDOS) of the 
zigzag GNR in the energy interval where only the lowest 
propagating state is occupied, see Fig. O A rectangular 
MB strongly modifies the lowest propagating state in the 
transverse direction. The wave function patterns in the 
barrier region can be easily understood from analysis of 
the corresponding patterns of Bloch states in the homo- 
geneous wire. The latter are shown in the right column 
of Fig. [SJ The lowest propagating state in the absence 
of the MB (Fig. [5] (a)) extends across the whole GNR 
at this energy. At B = 8 T, its probability density has 
a node at about 7nm away from the edge, and a local 
maximum is formed close to the center of the GNR. As B 
increases, this structure is pushed towards the edge of the 
GNR while its shape persists. A comparison of these pat- 
terns demonstrates that the wave function in the barrier 
region is directly related to the corresponding cigenstate 
of the homogeneous channel. Note that due to reflection 
on the barrier boundaries the edge state circulates inside 
the barrier region such that in this region |?/>| 2 has simi- 
lar amplitude near the upper and lower edges of the wire. 
We note further that the rapid oscillations corresponding 
to the wave functions on the sites belonging to the A and 
B sublattices are averaged out in the grey-scale plots to 
the left. 

The presence of this lowest propagating state appar- 
ently hampers the control of the carrier confinement in 
GNRs by MBs. It would thus be important to find a 
way to localize the lowest propagating state in gapless 
GNRs. There has been a lot of theoretical effort to ex- 
plore a way to open a bandgap in metallic GNRs. such 
as application of uniaxial strain, boron doping, and in- 



troduction of a line of impurities- 5 . To the best of our 
knowledge, no metallic behavior of GNRs with widths 
as studied here has been observed experimentally^. It 
was pointed out that the major discrepancies between 
experiments and theory may arise from the assumptions 
of perfect GNRs with a well-defined type of edge used 
in most theoretical studies. Experimental observations 
reveal that edge disorder is very significant on natu- 
ral graphite edges and etched GNRs. Theoretical stud- 
ies have shown that the edge disorder dramatically af- 
fects the transport properties and may turn the metallic 
ribbons into semiconductors^. Since the edge disorder 
is usually present in realistic GNRs, we investigate the 
transport properties in such GNRs subjected to a MB. 

In Fig. [HI we show the conductance of a zigzag GNR 
with edge disorder with and without a MB as a func- 
tion of Fermi energy The edge defects arc implemented 
by randomly removing 30 % of the atoms at the edges 
on both sides of the GNR, both inside and outside the 
magnetic barrier region. We first look at the conductance 
behavior in low energy regime. The characteristic feature 
is the appearance of conductance dips at the specific val- 
ues of the Fermi energy. Similar dips are also found in 
GNRs with additional bonds attached to the edges^. As 
the concentration of edge defects increases, the dips be- 
come more prominent and more zero-conductance dips 
appear, and their position and structure changes as the 
defect configuration is varied (not shown) . When the MB 
is activated, the position of some conductance dips, pre- 
sumably those originating from defects underneath the 
MB, move in energy while their amplitude is suppressed. 
The effect of the magnetic field is therefore to delocalize 
the lowest propagating state which have been localized 
by the edge defects. These results suggest that a mag- 
netic barrier can in fact be used to switch the conduc- 
tance in a GNR, but the mechanism differs from that one 
to be expected for magnetic barriers in two-dimensional 
graphene. Activation of the MB is able to delocalize the 
lowest propagating state in GNRs with disorder, thereby 
switching the conductance from zero to one. 



V. SUMMARY AND CONCLUSIONS 

We have provided an extensive theoretical study of the 
transport through graphene nanoribbons under the influ- 
ence of magnetic barriers. The magnetic barrier modifies 
the energy spectrum of the nanoribbon locally, which re- 
sults in an energy shift of the conductance steps. In ad- 
dition, multiple reflections along the transport direction 
between the entrance and the exit of the magnetic bar- 
rier generate Fabry — Perot resonances, the magnitude 
of which depends on the gradient of the magnetic field. 
These Fabry — Perot resonances are strongly suppressed 
in the case of magnetic barriers with smooth confinement. 
The lowest propagating state present in zigzag and metal- 
lic armchair GNRs is only weakly modified by magnetic 
barriers of realistic strengths. However, localization of 
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FIG. 5: Spatial distributions of local density of states in the zigzag GNR with N — 151 (left column) without a magnetic barrier 
(a), with a rectangular barrier of strength B — 8T (b) and B — 100 T (c) extending from x = — 60 nm to x = 60 nm. The gray 
scales (indicated by the bars on top) give the LDOS in units of t^nm . The Fermi energy is E F = O.Olto corresponding to 
the case when only the lowest propagating state is occupied. The corresponding probability densities of the lowest Bloch state 
(in arbitrary units) in an infinite homogeneous wire are shown in the right column. 
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FIG. 6: The Fermi energy dependence of the conductance of 
defective graphene ribbons with and without the rectangular 
MB of strengths 8T. The edge defect concentration is 30%. 



the lowest propagating state by disorder can be lifted by 
a perpendicular magnetic field, which offers a concept 
for magnetic barrier induced conductance switching in 
GNRs with disordered edges. 

In this paper we also present a novel method based 



on the Greens function technique for the calculation of 
the magnetosubband structure, Bloch states and mag- 
netoconductance of the graphene nanoribbons in a per- 
pendicular magnetic field. The non-trivial part of the 
method is the calculation of the surface Greens func- 
tion r, which typically requires very time-consuming 
self-consistent calculations. We, however, introduced 
a novel way to calculate the surface Greens function 
that does not require self-consistent calculations r^and 
where T is simply obtained from the solutions of the 
eigenequation of the dimension 2N x 2N (with N being 
the width of the nanoribbon). Utilization of this method 
obviously greatly facilitates computations, making the 
present method by far more efficient in comparison to 
the existing methods based on the self-consistent calcu- 
lations of r. The programming codes are freely available 
in the AIP EPAPS electronic depository^ 
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APPENDIX A: HOPPING MATRIXES U 

In this appendix we provide explicit expressions for 
hopping matrixes Un> , Eq. (|11[) . for armchair and zigzag 
ribbons in the Landau gauge A = (—By,0). The num- 
bering of slices and sites, r = within a unit cell 



J 



is given in Fig. ([2]), and the definition of the phases 
^r.r+A = 27r</v r+ A/0o an d the corresponding line inte- 
grals r , r +A are given by Eq. (j4|. In the expressions 
given below yi j stands for the y-coordinate of the site 



a. Armchair graphene ribbon 



where 0o,i;ij = —Byoja; 



02,1 = -k 



(Ui,o) iit = -toexp(i0o,j;i,j)6j,j>; 0o,i = 0+, 



/gi^l,l;2,l g^l,2;2,l 



gi#l,JV — 2;2,AT — 1 

gS#l,JV-l;2,JV-l g«#l,iV-l;2,JV 



V 



/ 



where <f>i j-2 j 
where 02,j;3,j 



-By 2 ja; 



(U 3 ,2) j;jl = -toexp(ie2,j;3,j)Sj,j>; 02,3 = 3 + 2 



/ gjS3,l;4 



04,3 = -to 



"3,1;4,2 gW3,2i4,2 



^3,Af-l;4,iV-l 
3 *^3,JV-l;4,iV P *03,iV;4,iV 



01,2 = 02+! 



J 03,4 = 



+ 
4.3 



where ^3 ,-.4 



f (»3j-a- 



^a 2 



1 = — -f ( J/3 ,a + -^a 2 ) ; and, because of periodicity, 



7 <ra,j;4j + i — "2 I ysj" ' 4 

05,4 = 01,0, 04,5 = 00.1 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



b. Zigzag graphene ribbon 



( u i,o)jj> = - i oexp(i6 l 0j . lji )(5 ij '; U ,i = 0+ o , 



where for odd j: 4>o,j-x,j = ~^B (yo.ja, + |a 2 ) , and for even j: (/>o,j;i,j = — ^rB (yo,jd — \a 2 ) 



V3 , 



where for odd j: 4>ij- t2 ,j 
periodicity, 



v3 
2 



(A6) 



(02,i)^v =-toexp(^ ljJ -. 2>i )^; 0i,2 = 2 + i, (A7) 
B (yija — |a 2 ) , and for even j: 4>ij- 2 ,j = ~ 2 ^ (yij' a + I fl2 ) > anc ^' because of 



03.2 ~ 01,0, 02,3 — 00.1- 
I 



(A8) 
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